Superconductors without an inversion center of symmetry: The s-wave state. 
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In materials without an inversion center of symmetry the spin degeneracy of the conducting band 
is lifted by an antisymmetric spin orbit coupling (ASOC). Under such circumstances, spin and parity 
cannot be separately used to classify the Cooper pairing states. Consequently, the superconducting 
order parameter is generally a mixture of spin singlet and triplet pairing states. In this paper we 
investigate the structure of the order parameter and its response to disorder for the most symmetric 
pairing state (A\). Using the example of the heavy Fermion superconductor CePtaSi, we determine 
characteristic properties of the superconducting instability as a function of (non-magnetic) impurity 
concentrations. Moreover, we explore the possibility of the presence of accidental line nodes in 
the quasiparticle gap. Such nodes would be essential to explain recent low-temperature data of 
thermodynamic quantities such as the NMR-Tj - , London penetration depth, and heat conductance. 



I. INTRODUCTION 

Early studies of superconductivity in materials with- 
out inversion symmetry addressed two-dimensional sys- 
tems, such as thin films, interfaces, and surfacesi2i2ii. 
The recent discovery of the non-centrosymmetric heavy 
Fermion superconductor CePtsSi^ has drawn attention 
towards bulk materials. Particular interest arose with 
the observation that the upper critical field of CePtaSi ex- 
ceeds the paramagnetic limiting field considerably. The 
simplest interpretation of this finding is in terms of spin 
triplet pairing. However, this is in conflict with the com- 
mon believe that the absence of an inversion center pre- 
vents electrons from forming spin-triplet pairsi. Indeed, 
inversion together with time reversal are the key sym- 
metries required for Cooper pairing. In a time reversal 
invariant system, the lack of inversion symmetry is con- 
nected with the presence of an antisymmetric spin-orbit 
coupling (ASOC). This can be represented in the single- 
particle Hamiltonian by a term of the general form 

ag k ■ S (1) 

where the vector function g k is odd in k (g k = —g-k) 
and a denotes the coupling strength 6 .. The ASOC is in- 
deed detrimental for most spin triplet pairing states as 
noted by Anderson 7 . However, it was also found that 
triplet states whose <i-vector lies parallel to g k would 
nevertheless be stable. Such a spin triplet state has the 
full symmetry of the crystal point group, we will therefore 
call it the S'-triplet state from now on. 

The interpretation of the absence of paramagnetism in 
terms of spin triplet pairing is not unique. Paramagnetic 
limiting is also drastically reduced for spin singlet states 
in the presence of ASOC 6 i 8 i 9 i 10 i n i 12 . In fact, the ab- 
sence of the inversion symmetry leads to a breakdown of 
the strict classification into even-parity spin-singlet and 
odd-parity spin-triplet pairing; these states are mixed, re- 
sulting in a state containing both componentsi2*i Since 
all but one spin-triplet states are suppressed by strong 
enough ASOC (i.e. a 3> krjT c , which is usually the 
case), we concentrate here on the stable S'-triplet state. 



This "high-symmetry" state mixes with the "s-wave" 
spin-singlet state (which also has the full symmetry of 
the crystal point group), since both of them belong to 
the same trivial A\ representation of C± v , the generating 
point group in CePtaSi. We will call this combined phase 
the " s-wave state" . 

The presence of ASOC leads to a splitting of the elec- 
tron bands by lifting the spin degeneracy. Thus, the 
discussion of the superconductivity is in some sense a 
two-band problem in this case. Assuming that supercon- 
ductivity is restricted to a single band the basis func- 
tion for irreducible representations of point group C± v 
have already been determined by various groups^^iii. 
These studies found that that the quasiparticle gap for 
the (most symmetric) A\ state would have the form 
A oc fc 2 + ky + cfc 2 . which is nodeless in general. In con- 
trast to these earlier works, we will examine here the full 
two-band situation. We show how the symmetry prop- 
erties of the pairing interaction and the distribution of 
the density of states on the two Fermi surfaces can in- 
fluence the form of the pairing state and, in particular, 
can introduce (accidental) line nodes in the quasiparti- 
cle gap which are not dictated by symmetry. Since these 
nodes are accidental their position in momentum space 
will vary both with temperature and magnetic fields. It is 
also interesting to see how the s-wave state is affected by 
non-magnetic impurities. These properties could be im- 
portant to obtain information on the complex structure 
of the pairing state in this material. 



II. MODEL WITH ANTISYMMETRIC 
SPIN-ORBIT COUPLING 

The basic model used to describe the conduction elec- 
trons in crystals without an inversion center can be writ- 
ten as 

k,s,s' 
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where { c ka ) creates (annihilates) an electron with 
wave vector k and spin s, & = [cr x , cr y , a z ) is the vec- 
tor of Pauli matrices and <j is the unit matrbiii. The 
band energy ^ k — e k — fi is measured relative to the 
chemical potential fi. The antisymmetric spin-orbit cou- 
pling (ASOC) term ag k ■ a is different from zero only 
for crystals without an inversion center and can be de- 
rived microscopically by considering the relativistic cor- 
rections to the interaction of the electrons with the ionic 
potential2ii&. For qualitative studies, it is sufficient to 
deduce the structure of the g-vector from symmetry 
arguments" and to treat a as a parameter. We set 
(9%) ~ 1j where () denotes the average over the Fermi 
surface. The ASOC term lifts the spin degeneracy by 
generating two bands with different spin structure. The 
normal state Green's function, 

G (k,ioj n ) = G + (k,iuj n )a + (g k • &)G°_(k, iu n ), (3) 
with 



G° ± (k,iu n ) 



1 



1 



± 



£k- a \9k\ iuJ n- Cfc + ot\g k \ 



_ (4) 
and g k = g k /\g k \ (\d\ = vt)i can be diagonalized into 
the components corresponding to the two bands, using 
the unitary transformation 

U{k) = cos(8 k /2) — ism(9 k /2)(cos4> k a- y — sm ^k&x), 

(5) 

where g k — \g k \(sm 9 k cos <j> k , sin^ sin(/> k , cos 9^) defines 
the angles 9 k and 4> k . This allows us to express, 

G (k,kj n ) = GS(fc,iw n )ffi(fc)+G§(fc,iw n )CT 2 (fc), (6) 
through the two Green's functions, 

f 



ei, 2 (fc)' 



(7) 



where the quasi-particle bands £1,2 (fc) = £fc ± a\g k \ are 
split by the presence of the ASOC. The spin structure in 
the bands is described by the 2x2 matrix 



°lA k ) = UxAkWljk) = l/2(<7 + (g k ■ <t))a, 
4 M = Ux, 2 (k)UUk) = l/2(a - (g k ■ <x)) A , 



(8) 



which act as projection operators in spin space. Since 
the spins of quasiparticles in the two bands have opposite 
directions for a given momentum k, the matrices satisfy 
the conditions, 



ai(k)a 2 (k) 
aUk) 





<7l, 2 (fc) (9) 
and are reciprocal under the inversion operation, 

<5-i, 2 (-fc) =<r 2 ,i(fe). (10) 



Symmetry considerations lead to the following form 
for <7fc in CePtsSi in lowest order expansion in the wave 
vector k: g k oc (k y , — k x , 0). The resulting spin structure 
is visualized in Fig^ where we assume for simplicity that 
the original Fermi surface is spherical. 



plane k — k 



plane k 





FIG. 1: Fermi surfaces for g k cc (k y ,~ k x ,0) as in 
CePtaSi. The arrows show the structure of the quasi- 
particle spin. Only along the z-axis the spin degeneracy 
is preserved. 



The pairing interaction is generally given by 
T-ipair = 2 X! X! 1/s i s ^ s 2 ;i i( fe ' fe ') c L 1 c -fc S2 c -fe'4 c fe' s ' 1 ' 

(11) 

where V aiB3 y a > (k, k') is the pairing potential with the 
symmetry properties, 



si so ,s's' 



(fe,fe') 



sos-i ,SU s' 



i-k,k') 



= -K 1S2 , si4 (fe,-fc') = K;4, S2S1 (fe',fe)- (12) 

In the absence of the ASOC the system is inversion 
symmetric in fc-space and the interaction factorizes in 
an orbital and a spin part. The conventional s-wave su- 
perconducting state belongs to the trivial representation 
A\ g of the inversion symmetric point group G. When 
the ASOC is turned on, the point group is reduced to 
the non-inversion symmetric subgroup G", whose trivial 
representation A\ does not have a definite parity. In fact 
A\ is compatible with both odd and even representations 
of G, i.e. with both {Ai g ,T u }, where T u is an odd parity 
irreducible representation in G,for which g k is a basis 
stated. Hence, the general form of the pairing potential 
relevant for the realization of the superconducting state 
belonging to Ax of G' has to involve both basis functions. 
Explicitly written, this leads to the form 



{g k •r) SlS2 f s t , s , +T SlS2 {g k , ■ f)^J } ,(13) 



where f s 



(ia v ) SlS2 , and, f ai< 



(icra v ) SlS2 . To 



avoid ambiguity we set V > 0, and ej + ef + e z m = 1. 
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The first term of Eqn. i|13fl is diagonal in the con- 
ventional s-wave pairing channel, the second one for the 
5-triplet pairing channel. The last term describes the 
scattering of Cooper pairs between the two channels, 
which is a result of the absence of inversion symmetry. 
The microscopic origin of this last term is explained in 
the Appendix In particular we show the existence of 
Dzyaloshinskii-MoriyaiZii& type of interaction for both a 
weakly interacting Fermi liquid and for a Hubbard model 
near half filling. This kind of interaction gives a contri- 
bution to e m . This can give rise to a large interaction 
when the derivative in the density of states is large. This 
is the case close to a Van Hove singularity. 



III. IMPURITY SCATTERING 

Effects of disorder are described by potential scattering 
of the quasiparticlcs, which in real-space representation 
is given by 



Hi = u{r- r t )iPl(r)il; s {r)dr, (14) 



Hi m p ^ 7~ti 



where u(r) is the potential of a non-magnetic impurity, 
which we consider rather short-ranged such that s-wave 
scattering is dominant. We are interested in the disorder- 
averaged normal and anomalous Green's functions G and 
F, 



Gx^r 



r>) = -(T r {Mr)%(r')}) 
r') = (T T {Mr)Mr')}) 
(T T {4(r)i>t(r')}), 



(15) 



with r = (r, r) and where the bracket denotes the ther- 
mal average. 

For the impurity average we use the Born 
approximation 2 ^, neglecting the possibility of more 
than two scattering events at the same impurity, which 
is legitimate if the potential is small in comparison with 
the characteristic electron energy scale e F (^f- Fermi 
energy, or analogue to the band width). 

Formally, the impurity scattering enters the self-energy 
of the Greens function of the normal, £ Gj Fig- Ela) , and 
the anomalous type, tip, Fig. |2Ib). Their mathematical 
expressions read 



£<?(£cj n ) 
t, F (iui n ) 



dk' 

ttNq J (2tt) 3 
r f dk' 



irN J (2?r) 3 



G(k',iujn) 

F(k',iu> n ), (16) 



where T = Trni mp Nou 2 is the averaged scattering rate, 
N n = (Ni + N 2 )/2 and iVi, 2 are the densities of state 
(DOS) of the two bands at the Fermi level. We intro- 
duced the impurity concentration and the s-wave scat- 
tering potential u 2 . 
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FIG. 2: The self energies contribution due to the im- 
purities scattering in the Born approximation, of normal 
type Eg (a) and of anomalous type £_f (b). 



The Gor'kov equations with these self-energy contribu- 
tions are formally analogous to those obtained for system 
with an inversion symmetry 20 



iuj n ) - T, G (ioj n ) ) G(k, iu n ) 

+ (A(k) + t F {iuj n )^ pi(k,iu n ) = <7 

(do lT (-fe, -iw n ) + s£(-»w„)) Ft(fc, iu n ) 

-(tt(k) + £ F (iUn))G(k,iuj n ) = 0. 

(17) 

The two-band structure of the normal state is more con- 
veniently handled, if we use the Green's functions G\ t2 
and iq,2) as is easily obtained by the unitary transfor- 
mation JSJ, 



( G Q l)=U\k)GU{k) 



(18) 



The phase factor — is introduced to cancel the phase 
dependence of the anomalous Green's functions on the 
two bandsii 

The pairing interaction defined in (|13fl fixes the form 
of the order parameter to A(fe) = i{ip + d (g^ ■ a-)}a y . ip 
and d are complex order parameters which can be inter- 
preted as the magnitude of the s-wave spin-singlet and 
respectively protected spin-triplet components. The par- 
ticular form of A(fc) prevents the existence of inter-band 
terms in the Gor'kov equations 

({Gl^k,^)}- 1 - £ G (z^ n )) G lt2 (k,iu n ) 

+ (A 1 . 2 (k) + Y, F (iuj n ))Fl 2 (k,iuj n ) = 1 
({G? i2 (-fc,-iw n )}- 1 + E G (-iw n ))F 1 t (2 (fc,iw n ) 

- (A 12 (fc) + E F (iw n )) G h2 (k,iuj n ) = 0, 

(19) 

where in this case 

T f dk' 



£ G (iw n ) 



27riVo J (2tt) 



{Gi(k',iu n ) + G 2 {k\iuj n )} 
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T, F (iLOn) 



and, 



dk' 



2ttN J (2?r) 



{Fi(k',iu> n ) + F 2 (k',iui n )} 



A h2 (k) = (< t p±d\g k \). 



(20) 



(21) 



Thus the Gor'kov equations are diagonal in the band in- 
dex. Interband effects are contained in the gap function 
and the scattering self-energies. The scattering off an 
impurity does not change the spin of a quasiparticlc and 
with the impurity-average a certain translational sym- 
metry is restored, such that the two bands to not mix in 
this approximation. Interband effects occur only through 
virtual processes. 

Introducing the modified gap functions A 1 2 (fc, iuj n ) = 
Ai,2(fc) + Si?(zcj„) and frequencies iCo n — zcj n — S^io;™), 
the solution of the two-band Gor'kov equations is given 
by 



and, 



Gl )2 (fc, ILUn) 



F lj2 (k,iw n ) 



(S£ + |Ai, a | a +# a )- 



Ai, 



|Ai, 2 | 2 + e? 2 ) 



(22) 



(23) 



and Ai 2(fcjiw n ) and u) n have to be determined self- 
consistently. Going back to the spin basis we find 

G(k,ioj n ) = G 1 (k,iuj n )a 1 (k) + G 2 (k,iu! n )a 2 {k), 
F(k,iuj n ) = {Fi(k,iu} n )a-i(k) + F 2 (k,ioj n )a 2 (k)}ia y . 

(24) 



IV. CHARACTERIZATION OF THE 
SUPERCONDUCTING INSTABILITY 

The gap functions of the two bands are obtained solv- 
ing the self-consistent equation 

A a (k) = -k B T J ^J2 V k,k' F ^ k '^- ( 25 ) 
The corresponding pairing interaction is determined by 

7j<5 si,s'. 



x U <s (k')uJ s ,(-k')D s .Ak'), 



(26) 



with D(k') =diag(— exp(— i<j!> fc /), — exp(i0 fc /)). Using i(T3|) 
and keeping only the contributions to pairing in the A\ 
pairing channel, it follows 

%,k' = ^{[ e s + e t \g k \\g k ,\}a +{e s -e t \g k \\g k ,\]a x 

- em [\g k \ + \g k >\] 6 z -ie m [\g k \ - \g k ,\] By) . 

(27) 



We parameterize the distribution of the normal state 
density of states in the two bands by constants N\ — 
N (l + Sn) and N 2 = N (l - 6 N )- To lowest order ap- 
proximation one may consider Sn ~ Q^o(£ = 0)/-^o- 

The critical temperature and the basic structure of the 
gap function A = {ip, d) which characterizes the super- 
conducting instability, follow from the solution of the lin- 
earized form of the self-consistent system of Eqs. (IB It 
IB3|) . Introducing Eqn. (|B1|> in Eqn. (|B3|) . and using 
standard summation techniques, we find 



1 



N V 
with 



{fi(e c ,k B T)Q l Q + f 2 (T,k B T)Q l r } 



(28) 



/i(e c , k B T) 

/ 2 (r,fc s T) 



ln{e c /2irk B T) + ln(4 7 ) 



2 2nk B T 



and, 



-e t 



MlSfcl) 



€-s 



(29) 



Qr 



H-«^(lflfcl> a )( I ' 



(30) 



Here we have introduced the digamma function ^S(z) de- 
fined by ^(z) = d/dz\n(z\), and the Euler's constant 
In 7 = C w 0.577. For simplicity we have assumed the 
same cut-off energy e c for both bands. 

First we characterize the superconducting instability 
of the clean system, i.e. r = 0. Eqn. H28JI has a 
non-trivial solution, if at least one of the eigenvalues Ai 
of the matrix Q l Q is positive. In this case the critical 
temperature follows the standard BCS relation k B T c = 
2e c j/n exp(-l/(N V\')), where A' = max^A,) > 0. 
Furthermore the nucleating form of the gap function fol- 
lows from Q l A' = A' A'. 

In this context we introduce the term of dominant 
channel to denote the channel responsible for the super- 
conducting transition, and we call subdominant the other 
channel characterized by A" — minj(Aj) and Q l A" = 
A" A". 

Solving the eigenstate problem Q l Ai = AjAj, we find 
the sets 



2A' = 
A' cx 



(e s + e t ) + 2e m 5 + A 

-(e s - e t ) + A 
2(e m - e t S) 



and 



where 



2A" = -(e s +e t ) + 2e m 5-A, 

£ii f-{e s — e t )-A 
^ ^ 2(e m - e t S) 

A = \J (e s - e t ) 2 + 4(e m - e s 6)(e m - e t S), 



(31) 
(32) 

(33) 
(34) 
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and 5 = <5w(|gi.|) to simplify the notation. Our interest 
lies on the characterization of the instability for all possi- 
ble combinations of the three components of the pairing 
potential l|13[) . From Eqn. I|31[) we determine the condi- 
tions for the superconducting instability and the form of 
the nucleating pairing state given by Eq. 

In order to conveniently display a phase diagram of the 
three parameters we use the condition e 2 s + e\ + e 2 n = 1 
and represent these parameters by the angles of spherical 
coordinates (e s = s- wave spin-singlet, e t = S spin-triplet, 
and e m = mixing): 

e s = cos(^„)sin(0„), 
e t = sin(<k,)sin(0„), 
e m = cos(0 v ). (35) 

The character of the interaction as a function of the two 
spherical angles (9 v ,<j> v ) is shown in Fig. 



71 
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I 1 

-3ti/4 71/4 6 57t/4 

FIG. 4: The character of the eigenvalues A', Eqn. f l3H . 
and A" , Eqn. f 1331 . as a function of the two spherical an- 
gles (0 V , 4>v)- If A' < we don't have any superconducting 
instability (grey domain with "no SC"). The white do- 
main is characterized by A' > 0, and A" < 0. If A" > 0, 
the system is characterized by two critical temperature 
(grey domain with " A" > 0" ) 
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FIG. 3: The character of the three interaction com- 
ponents (e s = s-wave spin-singlet, e t = protected spin- 
triplet, and e m = parity violating) as a function of the 
two spherical angles (9 V , <j> v ). The abbreviations " e a attr. 
" and " e s rep." mean that the s-wave component of the 
interaction is attractive (e s < 0) or respectively repulsive 
(e s > 0). 



Obviously an attractive interaction either for spin- 
singlet or -triplet pairing or both would yield a super- 
conducting instability. If both e s and et are repulsive 
(> 0), the mixing term plays the decisive role. From Eq. 
(|31|l we see immediately that an instability is absent for 
small \e m \ (grey domain with "no SC" in Fig. [3} : 



\e m \ < y/e s e t where e s , e t > 0. 



(36) 



On the other hand, a large enough value of \e m \ can trig- 
ger superconductivity, even with e s ,et > 0. The inter- 
parity scattering can lower the pairing energy and gen- 
erate positive eigenvalue of the matrix Q l Q . This mech- 
anism is formally analog to superconductivity driven by 
interband Cooper pair scattering in a multi-band super- 
conductor. 

A second positive eigenvalue, A" > is possible, if both 
spin channels are attractive (e s ,et < 0) (grey domain 
with "e s and et attr." in Fig. |3J). From Eqn. (|33|l we see 



that the following condition has to be satisfied: 

|e m | < ^e s e t and e s , e t < 0. (37) 

The parity-mixing term tends to suppress the second in- 
stability. 

It remains only to determine the form of the 2- 
dimensional order parameter A = (ip,d) based on Eqn. 
(|32|l . which nucleates at the superconducting transition. 
There are four typical forms of the order parameter 
("0, c?) , which will help us to structure the following dis- 
cussion (A = (ip,d)): 



oc 



P j Spin-singlet, A f oc ( ^ ) Spin-triplet. 



A,„j x ( 1 j 1st band , A m2 oc ( _ 1 



2nd ba(3fi) 



These correspond to the pure spin-singlet pairing A s , the 
pure spin-triplet pairing A t and two extreme cases of the 
mixed pairings A m j , and A,„2 • 

Before we start to discuss the conditions favoring dif- 
ferent order parameter forms, we consider the corre- 
sponding gap topologies on the two bands for the four 
cases in Eq. (|38H . To be concrete we discuss now the 
case relevant for the CePt 3 Si, ie. oc (k y ,— k x ,0). The 
gap on the two Fermi surfaces is given by 



Ai, 2 (0) = (^±d^/3/2\ sin(0)|), 



(39) 



where 9 is the polar angle in fc-space relative the z-axis. 
Fig. displays the two bands gap function Aj.,2 for the 
three different order parameter A s , A t , A,„i of Eq. 13811 . 
The case A m 2 can be obtained from A,„i by exchanging 
Ai and A 2 . 

In the cases of the pure spin-singlet phase A s or the 
pure spin-triplet phase A t , the two gap functions Ai and 
A 2 have the same amplitude, with the difference that in 
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FIG. 5: The gap amplitudes Ai,2 of the two non degen- 
erated bands for three particular directions of the order 
parameter in the 2-dimensional plane, spanned by ip, d. 
The directions corresponding to the two dimensional or- 
der parameter A„, A t , A mi are given by l|380 



Moreover the condition in Eq. I|36|) has to be satisfied to 
guarantee a superconducting phase transition. The first 
two conditions (flTl) and l|4"2|) are necessary to compensate 
the mixing of the spin-triplet with the spin-singlet pairing 
channels induced when the DOS are different for the two 
Fermi surfaces <5 7^ 0. Furthermore taking into account 
of the condition i|3(j|) we see that surprisingly A s and A 4 
can also nucleate in presence of repulsive spin-singlet e s 
and spin-triplet e* potentials. 

If e s = e t the conditions l|43|l and H44|l indicate that 
the order parameter is either A m i or A r „2. For Sn = 
the sign of the parity-mixing potential e m determines in 
which of the two order parameter appears. 

% r n 



Idl > ! Idl < lyl 




the spin-triplet case point-nodes appear along the z-axis. 
For A m i the gap on the first Fermi surface is much big- 
ger than on the second, Ai > A2. Furthermore the gap 
function in the second band A2 has two line-nodes per- 
pendicular to the z-axis. The presence of line- nodes in 
one of the two non-degenerate bands is not a particular- 
ity of the mixed states A ml , but a general property of all 
states with |d|Ci > and \ip\ ^ 0, where G\ = \/S/2 
only in case of a spherical original Fermi surface. How- 
ever in general the topology of the nodes in the two gap 
functions belonging to the trivial representation A\ of 
the point group C± v is characterized by 



\d\ 



= 2 point-nodes 

< C\ 2 line-nodes for one band _L z-axis (40) 
> C\ no nodes 



Now we turn to the conditions under which the dif- 
ferent forms of the order parameter nucleate. We expect 
that the order parameter appears as pure spin-singlet A s 
or a pure spin-triplet A t , if the spin-singlet potential e s 
or the spin-triplet potential e t are attractive and domi- 
nant, respectively. This is the case if the density of states 
on the Fermi levels is equal Sn — and the parity- mixing 
interaction e m is absent. 

From Eq. (|32[1 the necessary conditions can be derived 



A = 


A s 


for 




= e t S, e s < e t , 


(41) 


A = 


A, 


for 


&m 


= e s S, e t < e s , 


(42) 


A = 


A ml 


for 


e s 


= e t , e rn > e t S, 


(43) 


A = 


A m2 


for 


e s 


= e t , e m < e t .5 


(44) 
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FIG. 6: The direction of the nucleated order parameter 
A = (ip,(T). The grey lines show different possible lines 
for A = A s , and, A = At. The black ones are obtained 
for 5 = 0.3. 

Fig. El shows the phase diagram of the different order 
parameter forms resulting from Eq. (J3SJ. The differ- 
ent black lines are plotted choosing <5 = 0.3. The grey 
ones show how the lines corresponding to the pure spin- 
singlet A s and pure spin-triplet A t are moved varying 5. 
For general combinations of the three components e s , et, 
e m the angles 6 V , 4> v lie in one of the four basic domains 
delimited by the black lines shown in Fig. The order 
parameter nucleates with intermediary values of ip, and 
d. 



V. THE EFFECT OF THE DISORDER 

We now turn to the superconducting instability in a 
disordered system and analyze the behavior of the criti- 
cal temperature and the structure of the order parameter 
for different coupling parameters and impurity concen- 
trations. For this purpose we examine now Eq. 128|) 
for r > 0. First we would like to note that the ze- 
ros in the first column of the matrix Q l r are a conse- 
quence of Anderson's theoremSi, i.e. the conventional 
s-wave pairing state is not affected by non-magnetic im- 
purities. For more general states in our two-dimensional 
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order parameter space T c decreases with growing disor- 
der. We distinguish two basically different cases here. 
For the pure system either both eigenvalues A', A" are 
positive, or one of the two, A" is negative. In the first 
case there is a second lower (bare) critical temperature 
k B T^' = 1.14e c exp(-l/(iVoVA")), indicating the possi- 
bility of second superconducting phase transition. 



A. Case: A" = 

Before starting the discussion of this two general cases 
we consider the boundary situation with A" = 0. The 
interaction parameters e s , et, and, e m lie on a line given 
by, 



/e s e* with e s , e t < 0, 



(45) 



which correspond to the circle centered in <f> v = bir /A = 
— 37r/4 in the phase diagram of Figs. 0] and, HjJ The 
only eigenvector of the matrix Q r , with a non- vanishing 
eigenvalue corresponds to the form of A' given by in Eq. 
(|32[) . Hence in this case disorder would not alter the 
structure of the nucleating order parameter 



A' 



The instability equation 
f 



sign(e 
becomes 




(46) 



NnV 



= A7i(e c , k B T) + et(l - S 2 )f 2 (T, k B T) 7 (47) 



where et(l — S 2 ) is the eigenvalue of the matrix Q l r . We 
replace f x = hx{T c /T) + 1/N VX' and A' = -(e s + e t ) + 
2e m 5 in Eq. (|4T|> 



HTc/T) 



(e s + e t ) - 2e r , 
Using Eq. (|4(j|l we obtain 

1 



In 



where 



n 



i 



v 



:f 2 (T,k B T), 



«(7.*), 



(V; + Sd) 2 
d 2 {l-5 2 )' 



(48) 



(49) 



(50) 



with t = T/T c and 7 = T/k B T c . The effect of impurity 
scattering affects T c as in unconventional superconduc- 
tors in general 22 . However, a distinctive point is the pres- 
ence of the pre-factor 1/(1 + 77). A similar result was ob- 
tained for the s+g-superconductivity in borocar bides 23 . 

In order to visualize the behavior of the onset temper- 
ature of superconductivity depending on impurity con- 
centrations, we introduce the normalization with respect 
to the initial slope of T c -reduction. 
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{-dn imp /dt) t=1 (-dj/dt)^! 



(51) 



Fig. shows the evolution of t as a function of n' for 
different values of 77 when the subdominant channel has 
A" = 0. 




FIG. 7: The critical temperature as a function of the 
impurity concentration reduced so as to have a slope of 
— 1 at t = 1 for different value of rj. The value of 77, see 
eqn. 15UII . This result is valid only for A" = 0. 



For 77 > the superconducting instability extends to 
large impurity concentrations n' with the asymptotical 
behavior of the critical temperature 



0.88 rr'' -t 



4(1 + 77) 



(52) 



Hence we observe a variation of the robustness against 
non-magnetic impurities. For 77 = the standard behav- 
ior of an unconventional superconductor is found, while 
for rj — > 00 the non-sensitivity against disorder analogous 
to a conventional superconducting phase is realized. The 
latter case coincides with the suppression of the spin- 
triplet component in the pure phase. The exponent 77 
can be interpreted as the product of two ratios 



rj = 771772 with 771^ 



(jViAi(fc) + iV 2 A 2 (fc)) 
JV 1)2 <Ai(fc) - A 2 (fc)) ■ 



(53) 



where 771^ correspond to the ratio between the s-wave 
spin-singlet component (numerator) and the spin-triplet 
component lying in the first and respectively second band 
(denominator). Once (iViAi(fc) + N 2 A 2 (k)) / the ro- 
bustness of the superconducting state against disorder is 
introduced by the presence of a s-wave spin-singlet contri- 
bution. For other unconventional pairing states which do 
not belong to the trivial representation Ai , the impurities 
contribution to the self energy of anomalous type is zero 
and we recover exactly the result obtained by Larkin 22 . 



B. 



Case: A" / 



In this case the eigenvectors of the two matrices Q l 
and Q l r are different, such that we need to diagonalize 
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the matrix appearing in Eq. I|28(l . The highest eigenvalue 
and its eigenvector read 



2A = 2A7x + Cl - C 2 + yj4 + C3, 

III, \/ C 2+ c 3- c 4 \ 



(54) 
(55) 



0.5 



where (ip',d') are normalized + d' 2 — 1). Further- 
more, the parameters c, are defined as 



ci = (l-<5 2 ) 



e t (e t - e s - 2e m c>) + 2e 2 
A 



/ 2 . 



y ti=16 






r)=4 






11=1 






\ ri=0.25 




. V N 


r|=0.01 











c 2 = A fx - 



(1 - S 2 ) / : e? + e s [45(e m - e t 5) - e s 



A 



C3 = 



A(e 2 m ~e s e t ){e m -e t 8) 2 (5 2 -If 
A 2 



/I j 



c 4 = Afx + ct(l - <5 2 )/ 2 , 

C5 



v/4(e m - e t< 5) 2 + (e t - e s + A) 2 



(56) 



FIG. 8: The critical temperature as a function of the im- 
purity concentration and for different values of rj. Those 
corresponding to the dominant channel, solid black lines, 
have been reduced so as to have a slope of — 1 at t = 1. 
The shadow black lines show the evolution of the sub- 
dominant channel. This result is obtained supposing that 
the subdominant channel is attractive and t" = 0.2. 



The equation determining the critical temperature can 
be expressed as 



NnV 



= 2X'fx{e c , k B T) + cx-c 2 + Jc 2 + c 3 . (57) 



When we use again the relation fx = ln(T c /T) + l/A r V"A', 
we can cancel the 2/NqV term appearing at the left hand 
side of Eq. Q57JI and reach a new convenient representa- 
tion, 



The solid line shows the onset temperature of super- 
conductivity as a function of n' . The dashed line corre- 
spond to the bare transition temperature of subdominant 
instability, as obtained from the square root appearing in 
Eq. (E). 

The presence of an attractive subdominant channel 
supports the survival of superconductivity under non- 
magnetic impurity scattering. The lowest limit for the 
critical temperature is given by 







2 








+ C1 


- 2 


2V. n (I) 


+ c a 









C2 - C 3 = 0. 

(58) 

We assume now that the second channel is attractive 
A" > 0, so that it is possible to find a second solution 
for the instability temperature t" = T"/T c . We use the 
relation fx = ln(T c /T) + m(t")A"/(A" - A') to simplify 
Eqn. (5SJ. 

This leads to an equation of second order in ln(£) of 
the following form which allows us to determine both 
transition temperatures (t = T/T c or t = T" /T c ) 

In (J) (ln(0 - / 2 ( 7 ,t) - In(t)) = J_/ 2 ( 7> i)]n(t''), 

where now rj is given by the form of the order parameter 
in the clean system 



lim t = t 1+ " . 



(61) 



Using Eq. Ij57(l to simplify Eq. I|55|l we find that the form 
of the order parameter nucleating at the critical temper- 
ature t is given by 



■h(t) sf 3 (t) 
o 1 



d' 



where we introduced the function fs{t) defined as 

(l + r7)ln(i) 



hit) 



(1 + t?) In(t) - ln(t") ' 



(62) 



(63) 



+ Sd') 2 
d' 2 il-S 2 )' 



(60) 



and (ip',d') are connected to the pairing potential via 
Eqn. (32| . 

Fig. [S] shows the solutions of Eqn. (|59|l for different 
values of rj and for t 1 ^ = 0.2 in the pure case. 



The spin-singlet component of the order parameter in- 
creases with disorder. In the limit of dirty system the 
character of the order parameter is purely spin-singlet. 
Note that the subdominant instability is in any case sup- 
pressed by the disorder. The larger 77 the stronger the 
suppression. 

The behavior of the spin-singlet component ip of the 
normalized order parameter as a function of the critical 
temperature t and for different values of rj is shown in 
the Fig. [5] The order-parameter of the clean system is 
supposed to be in the spin-triplet channel. 
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T|=16 



0.5 








11=0.01 



FIG. 9: The spin-singlet component ip of the normalized 
order parameter as a function of the critical temperature 
t and for different values of rj. This result is obtained sup- 
posing that the subdominant channel is attractive with 
t" — 0.2 and that the order parameter of the clean system 
is spin-triplet . 



The case with a repulsive subdominant channel, A" < 
0, can be easily treated within the Eqs. I|59I63[1 substi- 
tuting ln(<") with l/NoV e ff, where 



NoVefj NqVX 1 N VX" 



(64) 



From Eqn. I|59|l we find that superconductivity disap- 
pears for a sufficiently high impurity concentrations. The 
critical impurity concentration n' c is given by, 



0.88 e 



^_ Tf _ = 4(l + r / ) 



(65) 



Furthermore the spin-singlet component of the nucleating 
order parameter decrease with disorder evolving towards 



-6d' 
d' 



(66) 



at the critical impurity concentration n' c . 

It has to be noted that for almost every combination of 
the pairing interaction A" ^ 0. However the realization 
of one of the two scenarios described in this chapter can 
be observed only if | A" | is very close or respectively bigger 
than A'. This means than in general we will observe a 
behavior very similar to that described by A" = 0. 



spin-orbit coupling. In addition to the pure spin-singlet 
s-wave and spin-triplet pairing, the pairing interaction 
includes a parity mixing contribution corresponding to 
an inter-parity scattering of Cooper pairs. The combina- 
tion of the three types of pairing interaction (e s , et, e m ), 
and the distribution of the density of states on the two 
non-degenerate bands <5, determines the form of the order 
parameter, represented by a singlet and a triplet compo- 
nent, A = (?/>, d). For the case of CePtsSi, we found that 
if \d\ > \tp\ and ^ , then one Fermi surface would 
have accidental line nodes perpendicular to the z-axis in 
the quasiparticle gap. The position of these accidental 
line nodes will in general depend upon both temperature 
and magnetic fields. 

The two-component structure of the order parameter 
allows for two distinct pairing channels, a dominant and 
a subdominant gap. We have shown that depending on 
the properties characterizing the two channels, the bi- 
phase is affected by non-magnetic impurities in different 
ways. In all cases the dependence of T c on the impurity 
concentration is very characteristic and could be used to 
establish the realization of the ^i-phase and the prop- 
erty of the subdominant pairing state. We would also 
like to mention that the presence of line nodes may leave 
its specific low-temperature features in various thermo- 
dynamic quantities such as the specific heat, the London 
penetration depth and so on. Recent experimental data 
of the nuclear magnetic relaxation rate T^UiLSi and the 
London penetration depth give evidence for the possible 
presence of line nodes. Moreover, the presence of two in- 
stabilities in the linearized gap equation, dominant and 
subdominant, provide the possibility that two supercon- 
ducting phase transitions could appear. From our study 
we conclude that such additional transitions would be 
observed in rather pure samples only, since the subdomi- 
nant instability is very sensitive to disorder. This is com- 
patible with the fact that in the dirty system eventually 
only the conventional pairing component would survive, 
while all alternative pairing channels are suppressed. 
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VI. CONCLUSION 

In this paper we characterized the superconducting in- 
stability towards a state belonging to the trivial repre- 
sentation Ai in a material without an inversion center 
and strong antisymmetric spin-orbit coupling. The corre- 
sponding pairing state involves spin-singlet s-wave pair- 
ing as well as a spin-triplet component specific to the 



APPENDIX A: ORIGIN OF MIXED 
SINGLET-TRIPLET PAIRING INTERACTION e m 

In this Appendix we consider a microscopic origin for 
the mixed singlet-triplet pairing interaction e m given 
in Eq. 13. In particular, we show that that the 
Dzyaloshinskii-MoriyaiZJ£ (DM) magnetic interaction 
gives a contribution to e m . Such an interaction is well 
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known to exist for magnetic systems that break inversion 
symmetry. We then show how this interaction can arise 
from the single particle Hamiltonian given in Eq. 2 for 
both a weakly interacting Fermi liquid and for a Hubbard 
model near half filling. Note that the DM interaction is 
not the only microscopic interaction that contributes to 
e m . Nevertheless, this discussion suffices to show to e m 
is not zero in general and in some cases may provide a 
substantial contribution to the superconducting conden- 
sation energy. 



DM Interaction 



The DM interaction can be written 



where D(q) is a real vector that satisfies D(q) = 
—D(—q). Invariance of H B m under point group opera- 
tions leads to the constraint RD B q = Dq, where R is the 
proper part of the rotation (R — Det(R) x R). Note that 
g(k) satisfies the same symmetry relation (Rg R ^ — g k )- 
Consequently, the two vectors D (q) and gq arc not 
orthogonal. Extracting the pairing contribution from 
Eci. lAll leads to the following mixed singlet-triplet pairing 
interaction 

[-D(k k') ■ f ] siS2 fj, s , + r slS2 [D(k - k') ■ t]\ s , . 

(A2) 

If we assume the form D(k — k') — e m (g k > — g k ) and 
impose symmetry constraints that arise from Pauli ex- 
clusion, we arrive at the mixed parity term in Eq. 2. As 
will be shown below, this form for D can be justified. 



2. DM Interaction in a weakly interacting Fermi 
liquid 

The DM interaction is known to exist in materials 
without inversion symmetry. This implies that the DM 
interaction should arise as a consequence of the existence 
of g k in the single particle Hamiltonian. Here we calcu- 
late the contribution of g k to the DM interaction through 
a calculation of the spin susceptibility. In the normal 
state the static spin-susceptibility is given by 

,uj n )a J G(k + q,u! n ))} 

k 

(A3) 

Upon using the form for the normal state Green's func- 
tions and carrying out the trace over the spins it can be 
shown 

Xij(q) = -2n 2 B k B Tj2J2 { [M 1 - 9k ■ 9k+q) 



+ (e* • 9fc)(ej ■ g k+q ) + (e* ■ g k+q )(e 3 ■ g k ) 

+ ife x ej) ■ (g k+q - g k ) G 1 (k,iv n )G 1 (k + q,w n ) 

+ [M 1 - 9k ■ 9k+q) + (e* • fffc)(e,- ■ g k+q ) 

+ & ■ 9k+q)(^J ■ 9k) - x ^) ' {gk+q ~ fffc) 

■ G 2 (fe, uj n )G 2 (k + q, uj n ) 

+ [<%(! + 9k ■ 9k+q) - fa ■ 9k) fa ■ 9k+q) 

- fa ■ 9 k+q )fa ■ 9k) ~ l fa x ^j) ' (9k+q + 9k) 

■ Gi(fc, w„)G 2 (fc + q,oj n ) 

+ [M 1 +9k 9k+q) fa ' 9k) fa ■ gk+q) 

~ fa-9 k+ q)(e r 9 k ) + ifa x Zj) ■ (9k+q+9k) 

+ G 2 (fe, w„)Gi(fc + q,w n )} 

The DM interaction is antisymmetric under interchange 
of i and j : 

D(q) = -2^ B k B Tj2{[ G l( k ^n)Gi(k + q,U} n ) 
k,uj n 

- G 2 (fe,w„)G 2 (fe + q,u n )] [9k+q 9 k ] 
+ [G 2 (fc, w„)Gi(fc + q,uj n ) 

- Gi(fc,w„)G 2 (fc + q,u n )] [gk+q+ 9k\] ■ 

(A4) 

Carrying out the sum over Matsubara frequencies, ex- 
panding the expression for D(q) to linear order in q one 
finds 



D(q) = 2^ B J2 
k 



Vn 2 • q Vni • q 



gq 



V6 q V£i • q 

m - n 2 



(6 - 6) 



r(Va + V6) 



q-Wg 

Vni + Vn 2 



(A5) 



where = £j(fe), = (e"* 4 + 1) _1 is the Fermi distribu- 
tion function for band i, and V = V^,. To linear order 
in a, this gives 



_,, , 8ur,a d 2 n(e k ) 

D(q)-^j:^f9-v k9k . 

k k 



(A6) 



For example if g k = kjkp (valid for a material with 
point group O), this gives 



D(q) 



,r/ 9 



, 4 aJVi q , . , 

Kf o VVq Kf 



where Nq = ^^-|^ = o is the derivative of the density of 
states evaluated at the Fermi surface and %p the Pauli 
susceptibility. This is a perturbative form and the deriva- 
tive of density of states is considered at the Fermi surface 
in the absence of spin-orbit coupling. Estimates for a 
show that the value can be a considerable fraction of the 
band width. 
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3. DM interaction within the Hubbard for finite 
ASOC 

Here we derive the DM interaction in the strong cou- 
pling regime by finding the effective Hamiltonian that 
governs the low- energy excitations of the Hubbard model 
with ASOC in the large U/t regime. The technique used 
to extract the Hamiltonian for the low-energy excitations 
is similar to that used for the derivation of the t-J model 
starting from the Hubbard model. 

We choose as zeroth-order Hamiltonian the on site 
Coulomb repulsion 



u = uJ2 n * 



(A8) 



The eigenstates of U are Fock states in the Wannier repre- 
sentation. U divides the Fock space into two sub spaces: 



5' 
D 



[\n lh n Ul n 2h - ■ •> : Vriij. + n A < 1] 



[KT> n U> n 2T> ' 



3n 



2]. 



(A9) 



D contain at least one doubly occupied site, and S are 
all configurations with either one or zero electrons per 
site. The hopping 



T = T S +T a 



where 



ij,ss f 

T a = - iOLij -(T ss ,cl s C 3 s> 

ij,ss' 



(A10) 



(All) 



contains now an antisymmetric term (cty = — ctji) cor- 
responding to the ASOC contribution. This means that 
also the effective interaction resulting by the superex- 
change process has an antisymmetric part. This last cor- 
respond to the DM interaction 



H 



DM 



ij,s 1 —s 4 



u 



C \si C j$2 n j] n ji C \s 3 C iS4 



ij 



u 



■acij ■ (Si x Sj) . 



(A12) 



The interaction TC DM is a factor a/t smaller than the 
usual spin-spin interaction constant J. While a/t will be 
less than one, it is not necessarily small. 



APPENDIX B: SELF-CONSISTENT EQUATION 

The self consistent equations ()28(l requires the evalua- 
tion of uj n and A = (ip, d). The corresponding equations 
are given by the substitution of Eqs. i|22|) and (|23|l into 
Eqs. HU) and (J201, 

= w n +Q n YQ 1 (iuj n ), (Bl) 
ip + ipTQi(ioj n ) + dTQ 2 (iw„), (B2) 



1 



N V 



A = irkgT Q(^n)&, 



(B3) 



with 

Q(w n ) = 



Qi(iuJn) ( -e s 
1 - YQi(iuj n ) \ e m 

Q2{iu n ) f e m -e s 
1 - TQx(iw n ) V ~ e t e m 
TQl{iuu n ) 



1 - TQi(iu> n ) 



c rn 

-6t 



(B4) 



where we introduced the functions 
1 + S N " 



Ql{lLO n ) 



(1-S N ) 



\J&1 + H + d\g k \\ 2 , 


Ju% + \$-d\g k \\' 


0- + 5N)\g k \ 


(l-S N )\g k \ 


^n + \i> + d\g k \\ 2 A 


Jcjl + \^-d\g k \\> 


(l + ^)l9fe| 2 


(l-5 N )\g k \ 2 



^ 2 + \4, + d\g k \\i ^+W Z d\g k \\ 2 , k 

(B5) 



The parameter Sn fix the distribution of the DOS at 
the Fermi level of the two bands, N% — N (l + Sn) and 

^2=JV (1-^). 



1 V. M. Edelstein, Sov. Phys. JETP 68, 1244 (1989). 4 O. Dimitrova and M. Feigelman, Sov. Phys. JETP Lett. 

2 L. P. Gor'kov and E. I. Rashba, Phys. Rev. Lett. 87, 78, 637 (2003). 

037004 (2001). 5 E. Bauer et at, Phys. Rev. Lett. 92, 027003 (2004). 

3 V. Barzykin and L. P. Gor'kov, Phys. Rev. Lett. 89 227002 6 P. A. Frigeri, D.F. Agterberg, A. Koga, and M. Sigrist, 
(2002). Phys. Rev. Lett. 92, 097001 (2004). 



12 



7 P.W. Anderson, Phys. Rev. B 30, 4000 (1984). 

8 L. N. Bulaevskii, A. Guseinov, and A. Rusinov, Sov. Phys. 
JETP 44, 1243 (1975). 

9 P. A. Frigeri, D.F. Agterberg, and M. Sigrist, New. J. Phys. 
6, 115 (2004). 

K. V. Samokhin, Phys. Rev. Lett. 94, 027004 (2005) 

1 V. P. Mineev, Phys. Rev. B 71 012509 (2005) 

2 R.P. Kaur, D.F. Agterberg, and M. Sigrist, Phys. Rev. 
Lett. 94, 137002 (2005). 

3 K. V. Samokhin, E. S. Zijlstra, and S. K. Bose, Phys. Rev. 
B 69, 094514 (2004) [Erratum: Phys. Rev. B 70, 069902(E) 
(2004)]. 

4 I. A. Sergienko and S. H. Curnoe, Phys. Rev. B 70, 214510 
(2004). 

5 V. P. Mineev, Int. J. Mod. Phys. B 18, 2963 (2004). 

6 G. Dresselhaus, Phys. Rev. 100, 580 (1955). 



I. Dzyaloshinskii, J. Phys. Chem. Solids 4, 241(1958) 
T. Moriya, Phys. Rev. 120, 91 (1960) 
Note that the in some cases gh cannot be attributed to 
a single irreducible representation. In this case the repre- 
sentation r~ is composite containing more than one irre- 
ducible represenation. We avoid this case here, since our 
aim is to give a description of pairing in CePtsSi where Y~ 
is irreducible. 

A. A. Abrikosov, L. Gor'kov, and I. E. Dzyaloshnskii, 
Dover, New York (1975) 

P.W. Anderson, J. Phys. Chem. Solids 11, 26 (1959). 
A. Larkin, Sov. Phys. JETP Lett. 2, 130 (1965). 
Qingshan Yuan et al, Phys. Rev. 68, 174510 (2003). 
P. F. N. Hayashi et al, pp. |cond-m at/0504176| 
K. Samokhin, pp. cond-mat/0504342 



